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Abstract 

Efficient slow and fast light fiber devices based on narrow band optical para- 
metric amplification require a strict polarization control of the waves involved 
in the interaction. The use of high birefringence and spun fibers is studied 
theoretically, possible impairments evaluated, and design parameters deter- 
mined. 
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1. Introduction 

It has been demonstrated that narrowband optical parametric amplifica- 
tion (NBOPA) is a superb technique for inducing slow and fast light (SFL) in 
optical fibers p|. Large group delay tuning, over selectable, wide frequency 
bands, makes this technique a very promising candidate for many envis- 
aged applications. Record tunable SFL delays of communication digital data 
signals were achieved in dispersion shifted fibers (DSF) [lj. In particular, 
experimental demonstrations were performed for 10 Gbit/s and theoretical 
predictions for 40 Gbit / s digital signals have been given [j, 0] . 

Theoretical studies resulted into a good understanding of the NBOPA, 
SFL process and its intrinsic limitations [l], 0] under the key assumption of 
an ideal homogeneous, isotropic fiber. However, real fibers are not likely to 
be homogeneous nor isotropic, and NBOPA gain and delay are affected by 
the longitudinal variations of the fiber linear propagation parameters such 
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as the zero dispersion wavelength (ZDW) fluctuations [5| and the random 
birefrigence 0]. 

Fluctuations of the ZDW cause gain broadening and thus a decrease in 
the achieved delay. The main contribution to ZDW shifts comes from changes 
in the fiber core effective area occurring during the drawing process 0,0], so 
ZDW shift is not time varying. Eventually, the ZDW shift can be measured 
(1, 0] and the gain broadening can be highly mitigated by selecting uniform 
samples, and by increasing the pump power. 

The polarization sensitivity of optical parametric interactions is a much 
more intriguing issue 0, [h|: the parametric gain coefficient is maximum 
when the pump and the signal have the same state of polarization (SOP), 
while it vanishes for orthogonal SOPs. So, maintaining a strict control of the 
pump and signal relative SOPs along the fiber is a critical issue for attaining 
reliable SFL devices based on NBOPA. 

The most widespread fibers (i.e. telecommunication ones) are not isotropic, 
though their birefringence is very low (here, they will be referred to as low 
birefringence - LoBi - fibers). The residual birefringence stems from the 
manufacturing imperfections (e.g. small asymmetries of the fiber core cir- 
cular section) and from the fiber operating conditions (e.g. stresses, bend- 
ing, twisting and temperature changes) due to environmental conditions; all 
these factors eventually break the polarization degeneracy of the fundamen- 
tal mode, and the fiber becomes birefringent. The residual, stochastic, low 
birefringence causes a random phenomenology known as polarization mode 
dispersion (PMD) Qjj]. The parametric interaction is very sensitive to 
PMD, that modifies the pump, signal and idler SOPs in a random fashion 
along the fiber 0. Detailed analyses of the effects of PMD on NBOPA 
gain and SFL delay can be found in refs. [H], 15, 16]. The effect is very 



pronounced and harmful because: a) it increases with the signal-pump fre- 
quency detuning, which is very large for NBOPA [l7j]; b) differently from 
ZDW fluctuations, which are deterministic, PMD is a random, time varying 
phenomenon ll|, [12 . 

In this contribution, two special fiber types that can be used to control 
the polarization of waves interacting in a NBOPA, SFL device will be theo- 
retically and numerically studied. The paper aims at presenting fundamental 
design information for improving the NBOPA, SFL fiber devices beyond the 
present state of the art. High birefringence (HiBi) fibers, in which the SOP 



is maintained [18L are an obvious choice to mimic the isotropic ideal case, 



but not the only one. Significant reductions of random polarization effects in 
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Raman 19 1, Brillouin 20] and parametric (2ll . 22| amplifiers, have been re- 
cently predicted when unidirectionally spun (US) fibers are considered. Here, 
the weight of possible negative effects in HiBi fibers and the positive effects 
of unidirectional spinning for on NBOPA-SFL are quantified. 

The paper is organized as follows. In section [2] the equations describing 
the NBOPA with polarized fields, in different fiber types, are introduced. 
The main features of propagation in LoBi, HiBi and US fibers will be also 
recalled in this section. In section [3] the performance in LoBi fibers will be 
assessed. The study of NBOPA, SFL in HiBi and US fibers will be carried 
out in sections 0] and [51 respectively. Finally, conclusions will be drawn in 
section EJ 



2. Propagation model 



The aim of this section is to provide a unified model to describe the 
NBOPA propagation for different fiber types, e.g. LoBi, HiBi and US. More- 
over, the main features of the propagation of optical signals in such fiber 
types will be also reviewed. 

Let us define \ A p (z)), \A s (z,t)), \Ai(z,t)), the Jones vectors of the pump, 
signal and idler waves respectively. For detailed definitions of the ket |) 
symbol and the bracket operators |)(|, (|), one can refer to 12 . 

In the undepleted pump approximation, neglecting the nonlinear effects 
of the signal and idler on the pump and by considering a continuous wave 
pump, the equations governing the nonlinear interaction of the slowly varying 
envelopes for the pump, signal and idler read [14|, [l5|, [ly, |23|, |24| : 



d\A p ) 

dz 
d\A s ) 

dz 
d\Aj) 

dz 



[C P + S P (\A P ))]\A P ), 

[C s + X s (\A P )) + U s {\A P ))\ \A S ) + T s {\A P )) I A*) (1) 
[£i + *i (\A P )) + TZi (\A P ))\ \Ai) + Ti (\A P )) \A* S ). 



The operators jC p , s ,i account for the linear propagation properties, while 
S p , X s i , 7Z S) i, T s ^i are operators that depend on the pump wave \ A p ), and take 
into account the nonlinear effects that are relevant for each wave, i.e. self- 
phase modulation, cross-phase modulation, Raman scattering and four-wave 
mixing. The linear operators are described in this section, as the different 
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fiber type (LoBi, HiBi, US) properties are dictated by these operators. Non- 
linear operators present very complicated structures, that are described in 
the Appendix [HJ for the sake of completeness. 
The linear operators are defined as: 

d 1- _ 1 - _ d 
£h = -OLh + jPh ~ Pih-T^ ~ j^Pi^h) ■ o- + ~5 h ■ a—, h = p,s, i. (2) 

In eqs. [2J ct p , s ,i are the loss coefficients, (3 p , s ,i = P(^ P ,s,i) the mean 
wavenumbers at the optical angular frequencies of the pump, signal and idler 
u PtSt i = 27rco/Ap iSi j satisfying 2u p = u s + u>i. The nonlinear phase matching 
condition to be satisfied is: 2f3 p — /3 S — Pi = A/3 = /3 2p (w s — oj p ) 2 + /3 4p {oj s — 
uip) 4, /12 = — 2jP , where [3 np is the n-th derivative of /3(u>) with respect to the 
angular frequency, calculated at cu p [ill], and P — (^(0)1^(0)) is the input 
pump power. The parameters used in the following simulations are: ZDW, 
A = 1.5423 /im; X p = 1.53 ptm; (3 2p ~ (3 30 (cu p - lu ); /3 30 = 1-14 • 10" 40 s 3 /m; 
(3 Ap = -5 • 1(T 55 s 4 /m; P = 1 5 W. Then, the NBOPA phase matched 
wavelengths are: A s ~ 1.3927 /im, A« ~ 1.7096 /im. 

If the reference frame (z, t) used in eqs. [1] is travelling at the signal 
group velocity v g (u s ), one also gets f3\h = 0, for the signal (h = s) and 
Pih = l/vg(ws) — ^-/ v g( UJ i) fo r t ne idler (h = i). 

The effects of fiber birefringence are accounted for by the last two terms 
of the operator ([2]). In the first term the Stokes vector /3(z,u) is the local 
birefringence vector that describes the birefringence at each point within the 



fiber; a is the vector of the Pauli spin matrices [12j. In the second term 



5h = d/3/duj calculated at u>h [23 



Results for an ideal, isotropic fiber, can be obtained by setting (3 = 
Vz, Vo>. In real fibers, the properties of the birefringence vector, which 
change from one type to the other, are very important. For this reason, they 
need to be specified in detail. 

For LoBi (unspun) fibers j3(z,u) = /3 un (z,uj) = /?2, 0] is a random 
vector. Its evolution along z can be obtained by means of the so-called 
random modulus model (RMM) [26j], i.e. its components are generated by 
the following Langevin equations: 

-Jj = -p(3i + vrji, i = 1, 2 (3) 

where (i — 1, 2) are Gaussian stochastic variables of zero mean and vari- 
ance Vp = u 2 /(2p) and rji(z) for i — 1,2 are independent, Gaussian white 
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noises of zero mean and unitary variance. Note that (3^{z) = is set in the 
numerical solutions, as is commonly assumed 26j| and experimentally verified 



in most cases [27J. This condition means that the fiber does not exhibit any 
circular birefringence. 

In the RMM, PMD is actually described by two length scales: the beat 
length Lb = 2n / (y/2vp) , and the birefringence correlation length Lp = 1/p. 
The former depends on the frequency (L b (uj) = uiqL b {ojq) fuS) and describes 
the length scale of polarization changes, while the latter is frequency inde- 
pendent (Lp = 9 m in our simulations), and accounts for the length scale 
of birefrin genc e changes. Both lengths contribute to determine the PMD 
coefficient [26(, hereinafter defined as D = a/ (At 2 )/L, where (At 2 ) is the 
fiber mean square differential group delay (DGD), and L is the fiber length. 
The DGD is defined as the time delay between pulses, at the same carrier 
frequency, launched along the two principal states of polarization (PSPs) 



28l | in LoBi fibers, and along the birefringence axes in HiBi fibers. The 



PSPs are defined as those input SOPs whose corresponding output SOPs 



are frequency independent at first order 28]. The typical coefficient D, in 
LoBi fibers, ranges from 10~ 2 ps/Vkm, for low PMD fibers, to more than 
KT 1 ps/Vkm for high PMD ones. 

The second type of fibers, HiBi, are for example realized by inducing 
internal stresses (PANDA fibers) or by making an asymmetric core (elliptical 
core fibers). In both cases the large intrinsic guide birefringence, voluntarily 
introduced during the manufacturing process, dominates over random effects. 
Then, in HiBi fibers the linear birefringence vector (3 has a predominant, 
deterministic, linear contribution ftmBi whose modulus is related to the DGD 
At by /3hibJ = coAr/(AL). Input SOPs parallel to birefringence axes 

(±PniBi = ±/3ffiBi/ \i3~mBi ) are polarization eigenstates, and so they travel 
unchanged through the entire fiber length. For this reason HiBi fibers are also 
referred to as polarization maintaining (PM) ones. Though the deterministic 
birefringence is overwhelming, random mode coupling along the fiber still 
exists. Furthermore, input misalignment may also result in a nonvanishing 
cross-polarized orthogonal SOP. The tolerance of the SFL technique with 
respect to SOP misalignement has been investigated and will be discussedin 
the following. Moreover, for a more realistic simulation of the unwanted 
effects, a random component, described again by the RMM, has been added 
to the deterministic part of the birefringence j3 = fimBi + Pran [U| • The value 
of D has been chosen to yield the typical polarization cross-talk ratio (PXR) 
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29j of commercial HiBi fiber, always better than 20 dB. As for the effects 
of input misalignments, the delay was calculated numerically exploring all 
possible linear input SOPs. 

Finally, US fibers have been considered; fiber spinning is a manufacturing 
process routinely performed while drawing fibers from LoBi preforms. In par- 
ticular, periodic spinning functions, in which the fiber is turned alternatively 
clokwise and counterclockwise, are often applied to reduce external stresses 
and fiber PMD . Unidirectional spinning, besides reducing the DGD, has 



been predicted to enhance the SOP alignment of optical signals at different 



frequencies, in particular in nonlinear fiber amplifiers [19|, [20j, |21|, |22j. The 
birefringence vector of a spun fiber can be obtained from that of the unspun 
case through the transformation /3(z,uj) = R 3 [2(j)(z)}/3 un (z, uj) where R 3 is a 
Mueller matrix representing a rotation around the third axis in the Stokes 
space (^3). For an US fiber the angle of rotation is given by the constant spin 
function <f)(z) = 2irz/p, where p is called the spin pitch. It has been shown 
3l| that when p 2 <C L 2 B , the polarization properties of the US fiber can be 



effectively described by a simplified model (SM). In the SM, the fiber can be 
represented by an equivalent birefringence vector, with a random linear and 
a deterministic circular component: 

p eq (z) = (^6(4 v^6(*), ~x) T , (4) 

where are statistically independent Gaussian white noises and the mod- 
uli of the linear and circular components are given by: 

2L F (np) 2 AnL F fi 
^~ L 2 B [p 2 + (4rrL F ) 2 Y X ~ p ' W 

The results of the SM are recalled here because in section they will con- 
tribute to explain the mitigation of polarization effects . 

For LoBi and US fibers, several hundreds statistical realizations of the 
stochastic processes, and subsequent integrations of eqs. [[J have been realized 
to calculate the mean gain and time delay. 

To conclude this section, let us remark that when propagation in bire- 
fringent media is considered, the group velocity cannot be uniquely defined, 



as observed by Haus [32|. Exceptions are represented by the special cases 
in which the input SOPs coincide with the PSPs, in LoBi fibers, or birefrin- 
gence axes, in HiBi fibers. In those two different value is found for 
the group velocity for each PSP or axis. In all other conditions, the standard 
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formula AT g = d((3z)/duj = z/v g (u) for the group delay looses its physical 
significance. Hence, in the numerical integrations, the group delay AT p has 



been evaluated as the first moment of the pulse as a function of time [15l . 133 
(the input pulse was Gaussian, 70 ps FWHM) 



t (A s (z,t)\A s (z,t)) dt 

AT g (z) = J — (6) 

(A s (z,t)\A s (z,t)) dt 



In particular, the difference between the arrival time when the SFL pump is 
on and when it is off, (AT = AT° n — AT°^), is calculated according to this 
definition. 



3. Low birefringence fibers 



Let us first briefly review the effects of PMD on SFL propagation in 



NBOPA 14, 15, 16 



For LoBi fibers the loss of alignment between signal 
and pump SOPs causes a reduction of the mean gain, which in turn translates 
into a delay reduction. For small PMD coefficients the mean delay (AT) can 
be calculated with the ideal isotropic case formula [3, 16], where the gain is 
replaced by the mean gain, i.e.: 



(AT) 




(7) 



1-^tanh 



where Ld = —Pip/Pip an d the mean nonlinear length (L^l) is related to the 
mean gain (G) by: 

(8) 



cosh" 1 (y(G) 



(L NL ) 

For larger random birefringence, strong pulse distortion sets in, and the 
delay decreases faster than the gain. Though mean quantities follow the ideal 
relation, random birefringence causes a large uncertainty in the actual delay. 
This fact can be easily grasped from fig. [H where the NBOPA gain versus 
delay, for many realization of an unspun fiber, for a PMD coefficient D = 
0.05 ps/ Vkm, is shown and compared with the case (squares) for D = (ideal 
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isotropic fiber). In the real case, PMD shifts the phase- matching condition 



14 , 1 1 61 ] , and therefore the gain is reduced, while eq. [7] is still valid. As we 
said, the delay standard deviation, represented by vertical red bars, is very 
large; around 20 — 30 ps i.e. almost 50 percent of the mean achieved delay. 
Finally, pulses are affected by a severe distortion. The large fluctuations of 
the delay are probably the most detrimental impairment introduced by PMD. 
So, regardless of LoBi fibers being the most commonly used ones, using them 
for making stable NBOPA based SFL devices is highly hampered. 

4. High birefringence fibers 

For SFL in HiBi fibers, the key parameter is power splitting among polar- 
izations. In the ideal case, since birefringence axes are eigenpolarizations, no 
cros-polarization coupling occurs if signal and pump are co-polarized, and 
aligned along ±(3 HiBi at the fiber input. Under real conditions, however, 
some coupling always exists because of input misalignment, either of the sig- 
nal or of the pump, with respect to the birefringence axes, and/or because 
of small imperfections in the fiber. Therefore, the robustness of the SFL 
scheme against input misalignements has been tested. The main results are 
summarized in fig. [2l 

The delay was first calculated as a function of the polarization misalign- 
ment between the (linear) signal SOP and the fast fiber axis (+/3mBi), along 
which the pump is suposed to be launched (empty circles). As we see, the 
delay is essentially unaffected for almost all signal SOP misalignements; it 
decreases significantly only when the signal is launched almost orthogonally. 
The explanation of the small change in delay is that a large polarization 



dependent gain is generated in this case 15]; therefore, the signal polariza- 



tion is attracted, by a nonlinear polarization pulling effect similar to that 



described for Brillouin and Raman amplification [3J, |35j), towards the di- 
rection yielding maximum gain. In practice, the portion of the signal pulse 
polarized along the minimum gain direction (slow axis) is weakly amplified, 
so that the pulse center of mass actually coincides with that of the powerful 
pulse component on the fast axis. When the pump and signal SOPs are or- 
thogonal (i.e., 90 degrees misalignment) there is no gain, and consequently 
no SFL effect. Hence, the delay tends to coincide with that due to linear 
birefringence. This has been verified by propagating the signal pulse without 
the pump (solid line); as the misalignment between the signal and the fast 
axis increases, more signal power is launched on the slow axis, and delayed 
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because of the different group velocities. A mirror-like behavior is observed 
if we launch the pump on the slow axis (squares in fig. [2]). In this case the 
delay is the sum of the SFL induced delay and that due to the change in the 
signal propagation axis. 

The effects of the pump input SOP misalignment are shown in fig. [2] 
(empty diamonds). As more pump power is launched on the slow axis, or- 
thogonally to the signal, the gain, and consequently the delay, decrease; both 
tend to zero when the misalignement tends to 90 degrees. Remarkably, there 
are certain input pump SOP's for which the delay becomes negative. This 
stems from the fact that the input signal wavelength is left constant, at 
the phase matching condition obtained when the pump is launched on the 
fast axis. But as the effective pump power decreases because of the mis- 
alignement, the phase matching condition shifts, and consequently the signal 
spectrum is now in a frequency band where the delay is negative. This expla- 
nation has been verified by calculating the delay with the ideal exact formula 
(i.e. eq. 9 of [la]); the results are presented in fig. [2] by the dashed curve. 

We may conclude that the delay reduction due to signal and pump input 
misalignment with respect to fiber axis is negligible, if the angle is less than 
10 degrees, an easy condition to satisfy in practice. 

Finally, we considered the random coupling, bearing in mind that this is a 
very small effect in good-quality HiBi fibers. Typically the PXR, i.e. the ratio 
between power on orthogonal axes at the output of 100m of fiber, is better 
than 20 dB. To evaluate these effects of the residual random coupling on the 
SFL delay, a stochastic component, obtained through RMM, was added to 
the deterministic linear birefringence vector, so that the total vector becomes: 
(3 = (3 mBi + Pran- The value of D was chosen such that the probability 
of getting realizations with a PXR in excess of 20 dB was very low. The 
results of a set of statistical realizations of this random process show that 
the maximum spread in the time delay is less than 1 percent of the mean 
value, very close to the ideal value. We conclude that the effect of the residual 
random coupling on SFL delay is negligible in HiBi fibers. 



5. Unidirectionally spun fibers 

Last, the case of US is considered. Fig. [3] shows a statistical set of 
realizations of a fiber with the same value of D as in fig. [T], but spun with a 
pitch p = 2 m. We see that the delay standard deviation is greatly reduced, 
and gets down to 5 ps. An impressive reduction of the standard deviation is 
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obtained by further decreasing the pitch, as shown in fig. HJ where p = 0.5 m. 
Note that there is no longer a significant difference with respect to the case 
of an ideally isotropic fiber. 

This remarkable result can be explained by the fact that in US fibers, as 
the spin pitch decreases, the equivalent deterministic circular birefringence 
and the random linear birefringence decrease. This is shown in fig. where 
the strengths of the random linear and deterministic circular components of 
the equivalent vector are calculated from eqs. [51 It is then clear that US 
fibers behave similarly to ideal isotropic fibers, when the spin pitch is short 
enough. The final outcome is that the pump and signal SOPs remain almost 
parallel, for all possible SOPs launched at the fiber input. 

The enhanced parallelism is illustrated also by figs. E] and [7] , where 
the mean value of cos(9 PtS ) (9 P;S being the angle between the signal and 
pump SOPs in Stokes space, at the fiber output) is shown for the two cases 
(unspun and US), as a function of detuning from the signal carrier frequency. 
Note that the alignment increases as the spin pitch decreses; moreover, it 
increses with pump power. The latter effect is, once again, an indication of 
the nonline ar p olarization pulling effect that we mentioned in the previous 
section 34, 35 



6. Conclusions 

The problem of how to control the SOP of the signal and pump, for slow 
and fast light in narrow-band, Raman-assisted, optical parametric amplifi- 
cation, has been studied theoretically and numerically. It was shown that 
standard, unspun telecommunication fibers can exhibit delay fluctuations, 
caused by polarization mode dispersion, that are too large in order to yield 
reliable slow and fast light effects for practical applications. 

To mimic an ideal isotropic fiber, two options have been explored: high 
birefringence fibers, and unidirectionally spun fibers. 

For HiBi polarization maintaining fibers, the typical achievable polariza- 
tion crosstalk ratio is sufficiently good to guarantee that an isotropic-fiber 
like delay is obtained. The effects of practical signal input misalignments 
are also under control, thanks to the nonlinear polarization pulling effect, 
which attracts the signal SOP towards the direction yielding maximum gain. 
Furthermore, small input pump misalignents, to be expected in practice, are 
not affecting the delay to a significant level. 
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Figure 1: Delay vs. gain from the numerical solutions of eqs. [T](L=lkm) for Pq: a) 1 W; b) 
2 W; c) 3 W; d) 4 W. Squares, isotropic fiber (D = 0); dots, LoBi fiber (D = 0.05 ps/Vkm); 
solid line, eq. [7] Bars define the gain and delay standard deviation; their crossing point is 
the mean value. 
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Figure 2: Delay vs. input misalignment. Signal - fast axis misalignment: Solid curve, 
no pump; Circles, pump aligned with the fast axis; Squares, pump aligned with the slow 
axis. Pump - fast axis misalignment: Diamonds, signal aligned with the fast axis; Dashed 
curve, analytical result. 
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Figure 3: Same as fig. [TJ but for a US fiber, p = 2 
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Figure 4: Same as fig. HJ but for a US fiber, p = 0.5 m. 




Figure 5: Strength of the SM equivalent birefringence vector components, (linear, 
dashed line) and \ (circular, solid line), as functions of the spin pitch p, as from eqs. [5] 
(L B = 10 m, L F = 9 m). 
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Figure 6: Alignment factor cos(9 p . s ) between the output signal and pump SOPs in Stokes 
space, as a function of the frequency detuning from signal carrier frequency for a LoBi 
fiber; dotted, dashed and solid curves are for Pq = 1, 3, 5 W respectively. 
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Figure 7: Same as fig. O but for a US fiber, p = 2m, 
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Finally, unidirectionally spun fibers were considered. For fast spinning 
(0.5 m), they have been shown to behave essentially like ideal isotropic fibers. 
In fact, compared to unspun fibers, the polarization mean alignment is highly 
enhanced along the entire length of the fiber. An additional advantage of 
spun fibers is that they maintain aligned any state of polarization, so the 
only care they require is to align signal and pump at the fiber input along 
the same direction. 

From a practical viewpoint, high birefringence commercial fibers are read- 
ily available, though customization might be necessary for tayloring the dis- 
persion properties to match the typical wavelengths of pump sources. As for 
unidirectionally spun fibers, spin pitches down to a few mm are feasible with 
present technologies. Therefore, such fibers could indeed be the most suitable 
tools for implementing reliable, wideband, slow and fast light devices. 
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8. Appendix 

The nonlinear operators of eqs. [1] are defined accordingly to the following 
expressions: 

„NR 



T 



where 



with 
and 



J- 



X 



[2(a p \a p )i + \a;)(a;\} + 
+j [(xjuu(o) + x£u«(o)) (A»IA»>i + xJi22i(o)|a;)(a;|] (9) 
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(10) 



(11) 

(12) 

(13) 
(14) 
(15) 

(16) 

(17) 
(18) 
(19) 



with 

xSn(n) = x5ia(n) + x5aa(n) + xf 221 (^)- (20) 

and 

= cj p — = uj s — Up. (21) 

The expression of the nonresonant components of the nonlinearity coef- 
ficient Xkimn an d of the resonant components of the nonlinear Raman sus- 
ceptibility Xkimnity can De found in [24j]. The effective areas A p and A s i t n, 
A s2 A2 read 



Cg>_ 2 
</ P 4 > 



4p = ThT' (22) 



j _ \Jpl\J S,ll / OQ v 

sM1 " w ( } 

As2 ' i2 " <i|7^ ' (24) 

where here angle brackets stand for integrals over the trasversal modal profiles 
fj = fj( x >y) U = Vi s ii)- 111 the numerical simulation of this paper it has 
been assumed A p = A s i t n = A s2t i2- 
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